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An  Extension  of  Blade  Element  Momentum 
Theory  to  Incorporate  Nonlinear  Lift  and  Drag  Coefficients 


Dewey  H.  Hodges* 

U.S.  Army  Research  and  Technology  Laboratories  (A  VRADCOM ) 
Ames  Research  Center,  Moffett  Field,  California 


An  expression  for  the  steady  induced  Inflow  velocity  is  obtained  from  Made  dement  and  momentum  theory 
considerations.  The  lift  coefficient  is  allowed  to  be  a  quadratic  function  of  angle  of  attack  and  the  drag  coef¬ 
ficient,  incorporated  through  a  perturbation  procedure,  can  be  any  arbitrary  function  of  angle  of  attack.  A 
principal  range  of  interest  is  identified  in  which  the  function  Inflow  versus  Made  pitch  angle  is  single  valued.  The 
resalts  reduce  to  the  classic  formula  found  la  the  textbooks  when  the  lift  coefficient  is  linear  and  drag  Is  ignored. 


Introduction 

THE  expressions  for  induced  inflow  velocity  in  hover  and 
vertical  flight  based  on  a  combination  of  blade  element 
and  momentum  theories  found  in  textbooks  such  as  Ref.  1, 
lack  the  generality  of  dealing  with  nonlinear  airfoil  lift 
coefficient  c,  versus  angle  of  attack  a.  When  generalized  to 
incorporate  nonlinear  c,  the  analysis  should  also  include  the 
influence  of  drag,  which  may  change  the  inflow  by  as  much  as 
20%  in  the  stall  regime.  An  additional  refinement  that  is 
desirable  when  the  formula  is  to  be  used  in  aeroelastic  stability 
analyses,  is  to  ensure  that  a  usei.ii  estimate  for  inflow  is 
available  for  both  positive  and  negative  thrust  and  in  either 
ascending  or  slow  descending  vertical  flight.  Such  a  formula 
has  not  been  published  and  it  is  the  purpose  of  this  note  to 
present  one.  This  kind  of  formula  will  avoid  having  to  solve 
for  the  inflow  iteratively  when  the  lift  coefficient  is  nonlinear. 
It  should  prove  to  be  helpful  in  aeroelastic  stability  analyses, 
such  as  Ref.  2,  when  the  influence  of  the  stall  parameters  on 
the  linearized  aeroelastic  stability  is  analyzed. 

The  procedure  herein  is  to  combine  expressions  for  the 
thrust  from  blade  element  and  momentum  theories  with  the 
assumption  that  the  lift  coefficient  is  quadratic,  over  some 
range,  in  angle  of  attack.  The  influence  of  drag  will  be  in¬ 
corporated  quite  naturally  as  a  perturbation  of  the  vertical 
component  of  induced  inflow  velocity.  The  effects  of  wake 
rotation  induced  by  lift  and  drag  are  neglected. 

Formulation 

From  either  blade  element  or  momentum  theory  it  is 
necessary  to  know  the  induced  velocity  to  specify  the 
operating  condition  of  the  rotor.  Consider  first  the  thrust 
from  momentum  theory  where  we  define  the  induced  velocity 
v,  and  vertical  velocity  V  positive  when  the  flow  is  down 
through  the  rotor  and  the  thrust  per  unit  blade  length  d  T/dr  is 
positive  when  air  loads  produce  an  upward  force  on  the  rotor. 
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The  momentum  balance  yields 


dT 

dr 


=  4fprv,  If,  +  FI 


(1) 


where  r  is  the  distance  from  the  center  of  rotation,  being  a 
maximum  of  R  at  the  blade  tip.  It  is  helpful  to  non- 
dimensionalize  the  variables;  thus. 


r=f,  If,  +  FI 


(2) 


where 


f,= 


v, 

a r 


df/dr 

4*prWR3 


(3) 


and  the  absolute  value  signs  ensure  that  thrust  r  (positive  up) 
and  induced  inflow  velocity  ti,  (positive  down)  have  the  same 
sign,  as  they  must.  A  plot  of  r  versus  f,  for  F>0  (climb), 
F<  0  (descent),  and  F=0  (hover)  is  found  in  Fig.  I.  Notice 
that  for  a  given  thrust,  the  induced  inflow  may  be  single 
valued  or  assume  three  different  values.  The  various  regimes 
are  identified  in  the  legend  of  the  figure.  In  the  case  of 
descent,  for  example,  when  n,  is  less  than  -  F/2  (i.e., 
f,  +  F/2<0)  but  f,>0,  the  flight  condition  is  autorotation 
and  the  rotor  is  operating  in  the  windmill  brake  state.  For 
-  K/2<f,<  -  f',  momentum  concepts  do  not  apply,  and, 
except  for  small  V,  when  0(>  -f'  they  are  not  generally 
accurate.  Notice  for  zero  rate  of  climb  that  the  results  are 
antisymmetric  about  the  origin.  Gessow  and  Myers  (Ref.  I) 
do  not  define  the  signs  in  a  rigorous  manner  for  these  unusual 
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Fig.  1  Threat  r  versus  induced  tallow  velocity  v,  for  various  rates  of 
dimb  V  based  on  Eg.  (2)  with  P-  -0.05,  0,  0.05.  (All  quantities 
dimensionless). 


conditions.  These  conditions,  though  certainly  ususual,  may 
be  encountered  in  experimental  testing  or  in  the  course  of 
iterating  to  solve  for  the  equilibrium  operating  condition 
when  blade  elastic  deflections  are  present  (Ref.  3). 

Now  let  us  write  the  thrust  from  blade  element  theory, 
including  drag,  assuming  that  the  total  resultant  wind  velocity 
with  respect  to  the  blade  element  l/is  given  by 

U=Qr/cos+  +  negligible  terms 


dimensionless  variables  with  e  so  that 


V—tv 
v ,  =  «X 


(7) 


This  yields 


XIX+vl  =  rc(  — *(X+v)crf+0(e'?)  (8) 


where  we  have  let  t2  =bc/8r.  Terms  of  order  t2  and  higher 
are  assumed  to  be  negligible.  Note,  however,  that  this 
assumption  breaks  down  if  f  is  small . 

The  angle  of  attack  is  given  by 

«=*-*=*-, an-'(^i) 

=  — r  <r0-R-X)  +  OUJ)=6t  (9) 

r 


where  the  pitch  angle  8  is  scaled  as  *0  and  the  angle  of  attack  a 
is  scaled  as  6t.  Now  let  X  =  X0  +  X;t  so  that 


5=60+6i(=  y  (r/3— v— X0-X/t) 

The  lift  and  drag  coefficients  in  terms  of  6  are  now 

f/(8)  +S;t~nf-(8o) 

do 

Crf(S)  —cd(60)  +®/*'rr  (®o) 
do 


(10) 


(ID 


Substitution  of  Eqs.  (10)  and  (1 1)  into  Eq.  (8)  yields 
dr* 

*,($„)  — *(v+X„)crf(6#)  = 


where  wake  rotation  induced  by  lift  and  drag  forces  is  (X^+rX,)  Ip+Xj  +  X/tl  +0(e*)  (12) 

neglected.  The  thrust  per  unit  blade  length  is  thus 


dT 

dr 


bpc  0 2r2 

y^-(C(COS0-CrfSin«) 


The  lift  coefficient  can  be  of  no  higher  degree  than  quadratic 
if  we  expect  a  closed  form  solution  to  Eq.  (12)  without  the  t 
terms.  Thus,  let 


bpcU2r2 

=  ^T(c'“Crftan*)  (4) 

where  b  is  the  number  of  blades,  c  is  the  local  blade  airfoil 
chord  length  at  radial  location  r,  and  the  inflow  angle  4  is 
given  by 


tan$  = 


y+vt 

a  r 


V+Oj 


(5) 


c/=c/„  +ci,a  +  ci1a2 


(13) 


Here  c,0  could  contain  the  effects  of  steady  coning  as  well  as 
airfoil  camber.  Now,  the  equations  for  ««  and  *'  in  Eq.  (12) 
with  c,  from  Eq.  (13)  are 


where  r*r/R. 

It  is  now  a  simple  matter  to  equate  the  thrust  parameters 
from  the  two  theories  in  Eqs.  (1)  and  (4)  to  obtain  an  equation 
for  6t .  In  dimensionless  variables  one  obtains 

r*t>,  Id,  +  9\ 

•  £  (6) 

where  f*c/R.  We  now  assume  that  0,  and  9  are  of  O(c) 
where  <  is  in  some  sense  small  compared  to  unity  and  scale  the 


[/-ySgn(X+v)jxj+  [r+  ['/  W-")] 

xsgn(X  +  »)  |x„-sgn(x  +  v)  +  (r0-») 

+  j  (r0-*-)>]«0  (14) 

*'  >  _  ~ (»,'*~ha)crf(8d)sgn(X-f ») 

'  2X,  +  v+(r,+2f,«e)sfn(X  +  »)  '  ' 
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The  solution  of  Eq.  (14)  is  then 

\>=(—  y  -sgn(A+i.)[y +y(r0— *)] 

±[[y  -ySgn(A+»»)J  +sgn(A+x) lft0 
+  rt i0+t30{rfl-y)]^^  [/  — ySgn(A+»>)J  (16) 


Now,  v,  can  be  expressed  from  Eqs.  (15)  and  (16) 
vi  =  t\=t(\0  +  e\l)=vio+vil 


with 


O-o  + 


KV-olMc,'  +c,2d)±Z 
l-(oc,2lt) 


t>ac. 


07) 


08) 


-0(6, 0  +  P)cA<*o) 

0  =  - - - 

*'  ±2Z 


where 


+  r«c(j  +6c,2(f8- F)J 


and 


be 

o-  y -  sgn(0,  +  V) 


09) 


(20) 


(21) 


Straightforward  numerical  evaluation  of  v,  for  representative 
values  of  bc/8w,  8,  V,  c,0,  c,( ,  cl},  cd,  and  r  yields  the  simple 
result  that  the  sign  of  the  quantity  Z  must  be  taken  the  same  as 
the  sign  of  v,  +  V.  The  choice  of  this  sign  is  commensurate 
with  the  thrust  and  inflow  being  of  the  same  sign.  Note  that  r 
must  not  be  small  and,  if  oc,}  >0,  we  have  a  lower  limit  on  its 
value  in  Eq.  (18). 

With  2 = Zsgn  ( 6,  +  V)  we  have 

0i+f'=(%  +  n[(-?£yf^]  (22) 


and  the  drag  correction  term  (underlined)  is  normally  smaller 
than  unity  so  that  sgn(u,  +  V)  =sgn(t/(#  +  F).  It  should  be 
noted  that  the  drag  coefficient  can  assume  any  functional 
form  without  precluding  a  closed  form  solution.  Changes  in 
sign  for  vi()  +  F  occur  at  8=8„  ,  where 


ci,’ci2*0  (23a) 


ch*0.cl}=0  (23b) 


Ifc(;  =c,}  =0,  vig  +  Fhas the signofc(fl.  Whenc.  ,c.  xOthc 
principal  range  of  interest  for  8  is  when  8<6„2(  a>  0)  or 
0>0cr2  (°‘<0).  Here  t),  may  be  regarded  as  a  single-valued 
function  and  8„}  is  given  by 


Of  course,  Z1  must  be  positive  for  all  cases  as  well. 

In  summary,  sgn ( v,  +  V)  =sgn (8-8  )  with  6„  given  in 

Eq.  (23)  and  sgn (8„}  -8)  =sgn  a  with  8„  given  in  Eq.  (24). 
The  induced  inflow  velocity  may  be  calculated  from  Eq.  (22) 
with 


,  ■  f,__*V-o(Hcl,+cl,8)+Z 
l-(oc,2/r) 


(25) 


It  is  not  necessary  to  explain  here,  in  detail,  all  possible 
behavior  of  the  solutions  of  these  equations.  It  is  noted  that 
results  from  the  above  formulas  may  be  infinite  or  complex 
for  certain  values  of  the  parameters.  Depending  on  the  nature 
of  specific  cases,  however,  such  problems  can  usually  be 
circumvented  with  little  difficulty. 

As  is  evident,  the  present  results  reduce  to  the  classical 
formula  as  given  in  textbooks,  for  example  Gessow  and  Myers 
(Ref.  1),  when  c,  is  linear,  cd  is  ignored,  and  the  sign  of  2  is 
taken  positive.  The  choice  of  the  sign  of  2,  however,  is  not 
adequately  discussed  in  Ref.  1 .  One  possible  extension  of  the 
present  work  is  incorporation  of  the  wake  rotation  component 
induced  by  lift  and  drag. 
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